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Tunnel magnetoresistance and interfacial electronic state 
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We study the relation between tunnel magnetoresistance (TMR) and interfacial electronic states 
modified by magnetic impurities introduced at the interface of the ferromagnetic tunnel junctions, 
by making use of the periodic Anderson model and the linear response theory. It is indicated that 
the TMR ratio is strongly reduced depending on the position of the d-levels of impurities, based 
on reduction in the spin-dependent s-electron tunneling in the majority spin state. The results are 
compared with experimental results for Cr-dusted ferromagnetic tunnel junctions, and also with 
results for metallic multilayers for which similar reduction in giant magnetoresistance has been 
reported. 



I. INTRODUCTION 

The tunnel magnetoresistance (TMR) of ferromagnetic 
tunnel junctions (FTJs) has attracted much interest due 
to its potentials for technological applications.Era To re- 
alize the high MR ratio and low resistance required for 
technological applications, more information on the elec- 
tronic states and tunneling proces3-.in FTJs is needed. 
Thus far, a simple Julliere's modeflcl for the TMR ratio, 
given as 



M R ratio = 



2P l Pr 
1 + PlPr 



(1) 



where Pl(r) is the spin polarization of the left(right)- 
hand side ferromagnetic leads, has frequently been used 
to analyze the TMR ratio. Although this expression well 
accounts for the experimental results of the TMR ratio 
once we adopt the values of the spin polarization P ob- 
served, the meaning of P is not sufficiently clear. In or- 
der to clarify the relation between the realistic electronic 
state and the TMR ratio, first principles and realistic 
tight-binding calculations have been performed to give 
higher TMR ratios than those observedflu Also, as a re- 
sult of the roughness of the .-amorphous-like barrier, the 
effects of resonant tunneling.0 scattering by disorder, and 
spin-flip tunneling da a to m agnetic impurities have been 
intensively studied. BBtaJld 

In spite of these intensive theoretical studies, the role 
of the interface is not yet completely understood. Recent 
experiments, however, have indicated the importance of 
the interfacial electronic states: Moodera et ai.ta showed 
that a formation of quantum well states at the interface 
by Cu spacer decreased the TMR ratio in a oscillatory 
way. Teresa et aZ£3 showed that the sign of P depends on 
the combination of the atomic species of the ferromag^ 
netic leads and the barrier materials. LeClair et alJt3 
demonstrated that the Cr-dusted interface strongly re- 
duced the TMR ratio of Co/Al-O/Co FTJs, whereas the 
Cu-dusted interface showed a much smaller effect, sug- 
gesting a close relationship between the electronic density 



of states at interfaces and the TMR ratio of FTJs. As 
for the decay of the TMR ratio due to the quantum well 
state, Zhang and Levyc3 proposed that a loss of the co- 
herence due to non-ideal structure of the inserted spacer 
may give rise to an exponentical decay of the TMR ratio 
as a function of spacer thickness. The difference between 
Cr and Cu spacers, however, appears even for 0.1 nm 
thickness of the spacer, which may indicate that a dis- 
orderd configuration of atoms at the interface must be 
properly taken into consideration. 

In this paper, we focus our attention on electronic 
states of the magnetic impurities at the interfaces, which 
have partly been argued by LeClair et al.fj and will show 
that the observed phenomena for Cr-dusted FTJs can be 
explained by the concept that the s-electron tunneling, 
which is dominant in the tunneling process, is strongly af- 
fected by changes induced by magnetic impurities in the 
electronic states at the interfaces. Below we adopt the 
periodic Anderson Hamiltonian to model the impurity- 
doped Co/Al-O/Co type FTJs, calculate the tunnel con- 
ductance and TMR ratio by using the Kubo formula, 
and show that the TMR ratio of FTJs made of Co or 
permalloy is actually sensitive to the majority spin elec- 
tronic states induced by magnetic impurities at inter- 
faces. Strong reduction in the giant magnetoresistance 
(GMR) or inverse GMR induced by magnetic imrHj£itieS| 
in metallic multilayers has already been reportedJisO'ta 
Such reduction in GMR ratios has been interpreted in, 
terms of the change in the spin asymmetry parameter .La 
We calculate the effects of magnetic impurities on the 
GMR ratio in a generalized periodic Anderson model 
(s — d model), then compare the effects of magnetic im- 
purities between TMR and GMR . 



II. MODEL AND FORMALISM 

In order to describe TMR and the electronic states of 
magnetic impurities in the interfacial Co layers in Co/Al- 
O/Co type FTJs, we must take into account the following 
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two aspects. First, the positive spin polarization of the 
tunnel current and the reasonably large TMR ratio of 
Co/Al-O/Co type FTJs must be properly reproduced in 
the model. The simplest way to satisfy these conditions is 
to use s-electron tunneling that is spin-polarized by spin- 
dependent s — d mixing. Second, it is well known that 
transition metal impurities such as Cr and V in Co metal 
have magnetic moments antiparallel to the Co moment, 
and form the so-called virtijaLbound state (VBS) in the 
majority spin state of Coo'Ej Because the formation of 
VBS results from the mixing of the impurity eHevel with 
the s-band of Co metal, we inevitably include the s-band 
in the model as well as the <i-band in order to describe 
the ferromagnetism of Co leads. Thus, the periodic An- 
derson model can include the basic features necessary to 
elucidate the electronic states and the characteristics of 
TMR. 

The Hamiltonian of the periodic Anderson model is 
given by 

H = H S + H d + H sd . (2) 

H s is the Hamiltonian for conduction electrons in a tight- 
binding version given as 

h s = -t s J2 c U- + E w *4c*> ( 3 ) 

ija ia 

where t s is the nearest neighbor (n.n.) hopping integral 
for s-electrons being independent of materials, v s is the 
energy level of the s-electrons, and a indicates the spin 
(T or I). The second term of eq.(2) indicates the d-levels, 
and is given as 

ff d = 5>£4^- (4) 

icr 

The last term indicates the s — d mixing, as 

ff«i = -7^(i4 + ^). (5) 

itr 

We neglect d — d hopping because the tunneling is domi- 
nated by s — s hopping, and assume that, by taking 7 = 
in the barrier region, only s-electrons tunnel through the 
barrier. Because we deal with transition metal impuri- 
ties, vf a represents both the (i-levels v% of the ferromag- 
netic leads and those of the impurities Vi mpa introduced 
at the interface. 

The tunnel conductance at the zero-bias limit is calcit. 
lated in numerical simulations using the Kubo formula.EHI 
We adopt a simple cubic structure for a finite-sized sys- 
tem; the cross section of the system includes 12 x 12 sites 
and the barrier thickness is 4 in units of the lattice con- 
stant a . Thereafter, the parameter values are taken in 
units oft 5 , such that v s = 0.0, vl: = -2.5, vf = -2.0, and 
7 = 1.5 for ferromagnetic leads. The barrier potential is 
taken to be 9.0. As can be seen later, a rather large value 
of 7 is required to obtain an appreciable TMR ratio by 
the s — d mixing. 
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FIG. 1: Calculated results of tunnel conductance as a func- 
tion of the Fermi level ef for parallel (P) and antiparallel (AP) 
alignment of the magnetization of the ferromagnetic leads. In- 
set: TMR ratio calculated from the tunnel conductance as a 
function of ep. 



III. CALCULATED RESULTS 

We first calculate the tunnel conductance without im- 
purities, Tp and Tap, in parallel (P) and antiparallel 
(AP) alignment of the magnetizations of the ferromag- 
netic leads, respectively. The TMR ratio is defined as 

MR ratio = Tp ~ TaF . (6) 
1 p 

Figure 1 shows the calculated results of tunnel conduc- 
tance as a function of the Fermi energy €f . As a result of 
the hybridization gap, Tp is zero between = —2.75 and 

— 1.80 for | spin channel, and between €f = —2.25 and 

— 1.40 for I spin channel. Consequently, Tap = 0.0 for 
—2.75 < £f < —1.40. The inset of the figure shows the 
results of the TMR ratio as a function of ep. The high 
TMR ratio below ep = —1.4 is due to opening of the 
hybridization gap, but there is no current in AP align- 
ment in this energy region. Above £f = —1.0, there is 
almost no TMR effect, because the current is carried by 
spin-unpolarized s-electrons. The physically meaningful 
values of the TMR ratio appear in a rather narrow range 
of energy, (—1.3 < ep < 1.1) where the s-electrons of 
the minority spin states carry less current than the ma- 
jority s-electrons due to the stronger s — d mixing at 
ep. Because the range in which a meaningful TMR ra- 
tio appears becomes narrower with decreasing 7, we take 
ep = —1.2 with a rather large value of 7 = 1.5 hereafter, 
which brings about an MR ratio of ~ 0.25, which is 
not unreasonable when compared with the experimental 
values. 

Now we introduce impurities into a ferromagnetic 
metal-layer adjacent to the insulating barrier. Because 
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FIG. 2: Calculated results of tunnel conductance in the par- 
allel alignment as a function of the impurity potential V5 mp t 
for ep = —1.2 with a mp = 0.1. Inset: the number of up spin 
electrons ni mp | (solid curve) and s-component of the density 
of states (broken curve) of a single impurity on a ferromag- 
netic surface as a function of Vi mp f . 



the transition metal impurities of the Cr type in which 
we are interested may form magnetic moments antipar- 
allel to the magnetic moments of the ferromagnetic leads 
made, for example, of Co, there is a relation such that 
^impT > Mmpi- Because of the relation, the up spin state 
of the impurities forms VBS, while the down spin state 
is almost occupied. Therefore, in order to realize the sit- 
uation, we take Vimpi = —2.5, which is slightly below 
vf = —2.0, and treat V; mp f as a variable parameter. In 
order to gain insight into the electronic state of the im- 
purity, we calculated the density of states (DOS) of a 
single impurity on the metallic surface, which value is 
a reasonable approximation of the impurity DOS at the 
interface of FTJs. We show in the inset of Fig. 2 the 
dependence of the number of | spin electron n; mp | of the 
impurity as a function of Vj mp f . With increasing T^mpt > 
"■impT decreases, indicating that the spin polarization of 
the impurity becomes negative for larger values of Vi mp f 
because ni mp i ~ 0.79. Because the drop in nj mp | be- 
comes stronger when Vi mp -|- approaches ep = —1-2, we 
find that the impurity DOS or VBS is located near ep. 
Although we did not carry out self-consistent calculations 
for the impurity state, the magnetic moment of the im- 
purity is about 0.4 per ci-orbital for VimpT = —0.5 and 
Vimpi = —2.5, which gives an on-site Coulomb interac- 
tion of ~ 5 in units of t s , corresponding to ~leV for 3d 
transition metals when we take t s ~ leV. It should be 
noted, however, that the spin polarization of the ferro- 
magnetic leads is rather small; this is a shortcoming of 
the present model, which was adopted to reproduce both 
positive spin polarization and a reasonable TMR ratio. 

The tunnel conductance of FTJs including impurities 
with concentration c; mp is calculated by averaging over 



FIG. 3: TMR ratio as a function of Vi mpT for e F = -1.2. 
Solid and open circles are the results for Ci mp =0.1 and 0.2, 
respectively. 



ten samples for different impurity-configurations. The 
root-mean-square deviation of the conductance is almost 
one order of magnitude smaller than the conductance it- 
self. Figure 2 shows the tunnel conductance Tpj in the 
parallel alignment as a function of V; mp j. It can be seen 
that the conductance decreases pronouncedly when Vi mp -|- 
exceeds eF- This is because the s — d mixing decreases 
the s-component of the majority spin state near ep as 
the VBS approaches ep, and blocks the s-electron tun- 
neling. In other words, an anti-resonance of the s-states 
appears at the interfacial layer where the impurities are 
introduced. The broken curve in the inset of Fig. 2 shows 
the s-component of the DOS of an impurity on a surface 
of an electrode. It can be seen the anti-resonant state 
appears when Vi mp ^ ~ ep, although the correspondence 
between the Pi mp f -dependences of the DOS and Tp is not 
perfect, probably due to the opening of the hybridization 
gap in the periodic Anderson model. 

As a consequence of the reduction in the | spin con- 
ductance in P alignment, the TMR ratio dramatically re- 
duces as Vimpi approaches ep, as shown in Fig. 3, where 
the results for c; mp = 0.2 as well as those for Cj mp = 0.1 
are presented. The results indicate that the TMR ra- 
tio decreases when VBS appears near the Fermi level 
and weakens the s-electron tunneling. The reduction of 
the TMR ratio due to magnetic impurities, however, is 
quantitatively stronger than the experimental value. The 
might be due to the rather strong s — d mixing chosen in 
this model. 

The relative position of VBS of magnetic impuri- 
ties in|-4fai*sjtion metals has been intensively studied 
before ;E3tilE3 VBS for Cr impurities is close to the Fermi 
level, whereas that of V impurities is higher than the 
Fermi level. For Cu impurities, on the other hand, V; mp j 
may be nearly the same with vi = —2.5, and does not 
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form any VBS. Therefore, we may conclude that impu- 
rities corresponding to Cr and V may bring about dra- 
matic reduction in the TMR ratio whereas Cu impurities 
do not. The results of both the conductance and TMR 
ratio are thus in_.good agreement with the experimen- 
tal observations.li-3 Further, the present model showed no 
reduction in the conductance of the f spin state when 
magnetic impurities were placed on the Co leads several 
layers apart from the interface. Rather, it slightly in- 
creased near VJmpT ~ 0. This-.rcsult is not inconsistent 
with the experimental result S3 



IV. COMPARISON WITH CIP-GMR 

In order to compare the results obtained above with 
similar effects reported for GMR, we performed simple 
calculations of GMR in trilayers doped with magnetic 
impurities near the interface. As a simple model of 
Co/Cu/Co type trilayers, we adopted the s — d model, 
the Hamiltonian of which is given by adding the following 
Hamiltonian to eq. (2): 



(7) 



because (i-electrons may carry the current in this case. 
The n.n. hopping integral for rf-electrons is taken to be 
0.1. Now the d-band is broadened by the (i-electron hop- 
ping, and we take a smaller value of the s — d mixing; 
7 = 0.3. 

The lattice structure is assumed to be simple cubic, 
and the thickness of each layer is /6/12/6/ atomic lay- 
ers in units of the lattice constant. Instead of the site 
representation used for TMR case, we use a mixed repre- 
sentation (i,k||) to treat the semi-infinite systems where 
I is an index of the atomic layers and k|| = (k x , k y ) is a 
Fourier transform of the sites on each atomic layer. The 
Green's function G a (z) = {z — Ha-)' 1 with z = e + irj 
is expressed by a tridiagonal matrix where the diagonal 
elements consist of 2 x 2 matrices of 



^(k||) 



"7 



"7 



(8) 



and the off-diagonal elements are alos given by 2 x 2 
matrices including the hopping integrals of s- and d- 
electrons. The d- level vf a depends on the layer, and is 
expressed as vf a in eq. (8). As in the TMR case, we 

-2.0 when I belongs to the 



take vf+ 



-2.5 and vf. 



ferromagnetic (Co) layers, and vf 



-2.5 when I 



belongs to the nonmagnetic (Cu) layers. Because a set of 
parameter values different from those for the TMR case 
is used for the GMR case, the position of £f is different 
from that taken for the TMR calculation. 

To reproduce the GMR effect in our simple model, 
we introduce the interfacial spin-dependent, 
caused by an intermixing layer at interfaces.! 



Because the d-level of the minority spin state of Co atoms 
is above the Cu ci-band, Co atoms intermixed with a Cu 
layer will form VBS in the minority spin state of the Cu s- 
band, resulting in spin-dependent scattering to produce 
the GMR. Here, we introduce the following self-energy 
into the s-component of the Green's function: 



C7 



-HA' 



(9) 



which represents the effects of Co atoms dissolved into 
the Cu layer. Here c indicates the concentration of such 
Co atoms and A is the broadening of <f- levels due to 
the s — d mixing, which is taken to be A = 0.1, keep- 
ing in mind the width of VBS calculated in the first 
principlesH3 We further take r\ = 0.001 to reproduce the 
spin-independent resistivity. Here, we note that the con- 
centration c itself does not have sufficient physical mean- 
ing, but C7 2 controls the magnitude of the self-energy. 
The conductivity parallel to the layers can be calculated 
using the Kubo formula without vertex corrections, be- 
cause there is translational invariance along the planes. 
The GMR ratios thus calculated increase with increasing 
c and reach about 60% for c = 0.5. 

After successfully reproducing the GMR effect in our 
simple model, we introduce additional magnetic impuri- 
ties into the adjacent layers of the intermixed interface. 
For simplicity, we use the same expression as that in eq. 
(9) for the self-energy for the majority-spin state of the s- 
band caused by the additional magnetic impurities, while 
c and vi a are replaced with c- lmp and Vi mpcr , respectively. 
As in the calculation of TMR, we take VJ mp | = vf, and 
treat Vi mp ^ as a variable parameter. The GMR ratios 
calculated for c = c- lmp = 0.5 with £f = —2.0 are shown 
in Fig. 4 (open circles) as a function of the impurity 
potential Vimpt- It can be seen that an inverse GMR 
can be realized when Mmpt ~ £ f- Open squares indicate 
the results for a trilayer with 6/6/6 atomic layers with 
ep = —1.9, and the results show that the qualitative fea- 
tures are not affected by the thickness of the trilayer. The 
GMR ratio becomes minimal when V5 mp | ~ ep, which 
characteristic is also found in TMR systems. The reduc- 
tion, however, may be brought about by an increased 
scattering of s-electrons in the majority spin band due 
to s — d mixing with the impurity states in addition to 
a reduction in the s-component of the DOS. The effects 
of magnetic impurities such as Cr and V in the GMR 
case are almost the same as those obtained in the TMR 
case, though the origin of the reduction in TMR ratio is 
due to the blockade of s-electron tunneling caused by the 
anti-resonance in the majority spin state of s-band. 



V. CONCLUSION 

In conclusion, by a qualitative and systematic study 
of the dependence of the MR ratio on the type of im- 
purities at the interface, we showed that both TMR and 
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^impT 



FIG. 4: Calculated results of the GMR ratio as a function of 
the impurity potential Vi mp | for 6/12/6/ and 6/6/6 trilayers, 
respectively. 



GMR ratios are affected strongly when the magnetic im- 
purities form VBS in the majority spin state near the 
Fermi level. The results are in good agreement with the 
observed ones. Although the effects of magnetic impu- 
rities such as Cr and V are quite similar on TMR and 
GRM, the reduction in TMR ratios is due to a blockade 
of s-electron tunneling, whereas that in GMR is mainly 
due to an increase of s-electron scattering in the majority 
spin state. More realistic and self-consistent calculations, 
however, are desired to obtain a quantitative understand- 
ing between the MR ratio and the interfacial electronic 
states. 
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We study the relation between tunnel magnetoresistance (TMR) and interfacial electronic states modified 
by magnetic impurities introduced at the interface of the ferromagnetic tunnel junctions, by making use of 
the periodic Anderson model and the linear response theory. It is indicated that the TMR ratio is strongly 
reduced depending on the position of the d-levels of impurities, based on reduction in the spin-dependent .s- 
electron tunneling in the majority spin state. The results are compared with experimental results for Cr-dusted 
ferromagnetic tunnel junctions, and also with results for metallic multilayers for which similar reduction in giant 
magnetoresistance has been reported. 



I. INTRODUCTION 

The tunnel magnetoresistance (TMR) of ferromagnetic tun- 
nel junctions (FTJs) has attracted much interest due to its po- 
tentials for technological applications. 1,2 To realize the high 
MR ratio and low resistance required for technological appli- 
cations, more information on the electronic states and tunnel- 
ing process in FTJs is needed. Thus far, a simple Julliere's 
model 3 - 4 for the TMR ratio, given as 



MR ratio = 



2P l Pr 
1+PlPr' 



(1) 



where Pl(R) is the spin polarization of the left(right)-hand side 
ferromagnetic leads, has frequently been used to analyze the 
TMR ratio. Although this expression well accounts for the 
experimental results of the TMR ratio once we adopt the val- 
ues of the spin polarization P observed, the meaning of P 
is not sufficiently clear. In order to clarify the relation be- 
tween the realistic electronic state and the TMR ratio, first 
principles and realistic tight-binding calculations have been 
performed to give higher TMR ratios than those observed. 5,6 
Also, as a result of the roughness of the amorphous-like bar- 
rier, the effects of resonant tunneling, 7 scattering by disorder, 
and spin-flip tunneling due to magnetic impurities have been 
intensively studied. 8-11 

In spite of these intensive theoretical studies, the role of 
the interface is not yet completely understood. Recent exper- 
iments, however, have indicated the importance of the inter- 
facial electronic states: Moodera et al. n showed that a for- 
mation of quantum well states at the interface by Cu spacer 
decreased the TMR ratio in a oscillatory way. Teresa et al. n 
showed that the sign of P depends on the combination of the 
atomic species of the ferromagnetic leads and the barrier ma- 
terials. LeClair et al. u demonstrated that the Cr-dusted in- 
terface strongly reduced the TMR ratio of Co/Al-O/Co FTJs, 
whereas the Cu-dusted interface showed a much smaller ef- 
fect, suggesting a close relationship between the electronic 
density of states at interfaces and the TMR ratio of FTJs. As 
for the decay of the TMR ratio due to the quantum well state, 
Zhang and Levy 25 proposed that a loss of the coherence due to 
non-ideal structure of the inserted spacer may give rise to an 



exponentical decay of the TMR ratio as a function of spacer 
thickness. The difference between Cr and Cu spacers, how- 
ever, appears even for 0. 1 nm thickness of the spacer, which 
may indicate that a disorderd configuration of atoms at the in- 
terface must be properly taken into consideration. 

In this paper, we focus our attention on electronic states of 
the magnetic impurities at the interfaces, which have partly 
been argued by LeClair et al., u and will show that the ob- 
served phenomena for Cr-dusted FTJs can be explained by the 
concept that the s-electron tunneling, which is dominant in the 
tunneling process, is strongly affected by changes induced by 
magnetic impurities in the electronic states at the interfaces. 
Below we adopt the periodic Anderson Hamiltonian to model 
the impurity-doped Co/Al-O/Co type FTJs, calculate the tun- 
nel conductance and TMR ratio by using the Kubo formula, 
and show that the TMR ratio of FTJs made of Co or permal- 
loy is actually sensitive to the majority spin electronic states 
induced by magnetic impurities at interfaces. Strong reduc- 
tion in the giant magnetoresistance (GMR) or inverse GMR 
induced by magnetic impurities in metallic multilayers has al- 
ready been reported. 16-18 Such reduction in GMR ratios has 
been interpreted in terms of the change in the spin asymmetry 
parameter. 19 We calculate the effects of magnetic impurities 
on the GMR ratio in a generalized periodic Anderson model 
(s - d model), then compare the effects of magnetic impurities 
between TMR and GMR . 



II. MODEL AND FORMALISM 

In order to describe TMR and the electronic states of mag- 
netic impurities in the interfacial Co layers in Co/Al-O/Co 
type FTJs, we must take into account the following two as- 
pects. First, the positive spin polarization of the tunnel current 
and the reasonably large TMR ratio of Co/Al-O/Co type FTJs 
must be properly reproduced in the model. The simplest way 
to satisfy these conditions is to use s-electron tunneling that 
is spin-polarized by spin-dependent s - d mixing. Second, it 
is well known that transition metal impurities such as Cr and 
V in Co metal have magnetic moments antiparallel to the Co 
moment, and form the so-called virtual bound state (VBS) in 



the majority spin state of Co. 20 ' 21 Because the formation of 
VBS results from the mixing of the impurity J-level with the 
s-band of Co metal, we inevitably include the s-band in the 
model as well as the <f-band in order to describe the ferromag- 
netism of Co leads. Thus, the periodic Anderson model can 
include the basic features necessary to elucidate the electronic 
states and the characteristics of TMR. 

The Hamiltonian of the periodic Anderson model is given 

by 



H — H v + Hri + H. 



sd- 



(2) 



H s is the Hamiltonian for conduction electrons in a tight- 
binding version given as 



H s = -t s Y^ clc j(T + J] v s clc itT , 



(3) 



>J<r 



where f is the nearest neighbor (n.n.) hopping integral for 
^-electrons being independent of materials, u s is the energy 
level of the ^-electrons, and cr indicates the spin (| or X). The 
second term of eq.(2) indicates the cZ-levels, and is given as 



z'cr 

The last term indicates the s - d mixing, as 
H s <t = -y^{c] a d io - + h.c). 



(4) 



(5) 



We neglect d-d hopping because the tunneling is dominated 
by s - s hopping, and assume that, by taking y = in the 
barrier region, only s -electrons tunnel through the barrier. Be- 
cause we deal with transition metal impurities, v d icr represents 
both the J-levels v d a of the ferromagnetic leads and those of 
the impurities Vi mpix introduced at the interface. 

The tunnel conductance at the zero-bias limit is calculated 
in numerical simulations using the Kubo formula. 10 We adopt 
a simple cubic structure for a finite-sized system; the cross 
section of the system includes 12 x 12 sites and the barrier 
thickness is 4 in units of the lattice constant a . Thereafter, the 
parameter values are taken in units of f, such that v s = 0.0, 
\A = -2.5, iff = -2.0, and y = 1.5 for ferromagnetic leads. 
The barrier potential is taken to be 9.0. As can be seen later, 
a rather large value of y is required to obtain an appreciable 
TMR ratio by the s — d mixing. 



in. CALCULATED RESULTS 

We first calculate the tunnel conductance without impuri- 
ties, Tp and Tap, in parallel (P) and antiparallel (AP) align- 
ment of the magnetizations of the ferromagnetic leads, respec- 
tively. The TMR ratio is defined as 



MR ratio 



- AP 



(6) 
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Figure 1 shows the calculated results of tunnel conductance 
as a function of the Fermi energy e F . As a result of the hy- 
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FIG. 1: Calculated results of tunnel conductance as a function of 
the Fermi level e F for parallel (P) and antiparallel (AP) alignment 
of the magnetization of the ferromagnetic leads. Inset: TMR ratio 
calculated from the tunnel conductance as a function of e F . 



bridization gap, Tp is zero between ep = -2.75 and -1.80 for 
f spin channel, and between e F = -2.25 and -1.40 for J, spin 
channel. Consequently, T A p = 0.0 for -2.75 < e F < -1.40. 
The inset of the figure shows the results of the TMR ratio as 
a function of e F . The high TMR ratio below e F = -1.4 is 
due to opening of the hybridization gap, but there is no cur- 
rent in AP alignment in this energy region. Above e F = —1.0, 
there is almost no TMR effect, because the current is carried 
by spin-unpolarized s-electrons. The physically meaningful 
values of the TMR ratio appear in a rather narrow range of en- 
ergy, (—1.3 < e F < 1.1) where the s-electrons of the minority 
spin states carry less current than the majority s-electrons due 
to the stronger s-d mixing at e F . Because the range in which 
a meaningful TMR ratio appears becomes narrower with de- 
creasing y, we take e F = —1.2 with a rather large value of 
y = 1.5 hereafter, which brings about an MR ratio of ~ 0.25, 
which is not unreasonable when compared with the experi- 
mental values. 

Now we introduce impurities into a ferromagnetic metal- 
layer adjacent to the insulating barrier. Because the transition 
metal impurities of the Cr type in which we are interested may 
form magnetic moments antiparallel to the magnetic moments 
of the ferromagnetic leads made, for example, of Co, there is 
a relation such that Vi mp f > Vi mp j,. Because of the relation, 
the up spin state of the impurities forms VBS, while the down 
spin state is almost occupied. Therefore, in order to realize 
the situation, we take Vi mp j, = -2.5, which is slightly below 
v d = -2.0, and treat V; mp f as a variable parameter. In order to 
gain insight into the electronic state of the impurity, we cal- 
culated the density of states (DOS) of a single impurity on 
the metallic surface, which value is a reasonable approxima- 
tion of the impurity DOS at the interface of FTJs. We show 
in the inset of Fig. 2 the dependence of the number of f spin 
electron «i mp1 - of the impurity as a function of Vimpj. With 
increasing Vi mpT , «i mpT decreases, indicating that the spin po- 




FIG. 2: Calculated results of tunnel conductance in the parallel align- 
ment as a function of the impurity potential V imp -f for e F = -1.2 with 
t'imp = 0. 1 . Inset: the number of up spin electrons n impT (solid curve) 
and s-component of the density of states (broken curve) of a single 
impurity on a ferromagnetic surface as a function of V impT . 



larization of the impurity becomes negative for larger values 
of Vimpt because n; mp | ~ 0.79. Because the drop in n; m pt be- 
comes stronger when V m pt approaches e F — - 1 .2, we find that 
the impurity DOS or VBS is located near e F . Although we 
did not carry out self-consistent calculations for the impurity 
state, the magnetic moment of the impurity is about 0.4 per 
(i-orbital for V m pt = -0.5 and V; mp j, = -2.5, which gives an 
on-site Coulomb interaction of ~ 5 in units of t 5 , correspond- 
ing to ~leV for 3d transition metals when we take t s ~ leV. It 
should be noted, however, that the spin polarization of the fer- 
romagnetic leads is rather small; this is a shortcoming of the 
present model, which was adopted to reproduce both positive 
spin polarization and a reasonable TMR ratio. 

The tunnel conductance of FTJs including impurities with 
concentration a mp is calculated by averaging over ten samples 
for different impurity-configurations. The root-mean-square 
deviation of the conductance is almost one order of magni- 
tude smaller than the conductance itself. Figure 2 shows the 
tunnel conductance r P1 - in the parallel alignment as a function 
of Vimpt- I* can b e seen mat tne conductance decreases pro- 
nouncedly when Vi mp f exceeds e F . This is because the s - d 
mixing decreases the ^-component of the majority spin state 
near e F as the VBS approaches e F , and blocks the s-electron 
tunneling. In other words, an anti-resonance of the s-states 
appears at the interfacial layer where the impurities are intro- 
duced. The broken curve in the inset of Fig. 2 shows the 
s-component of the DOS of an impurity on a surface of an 
electrode. It can be seen the anti-resonant state appears when 
Vmpt ~ f F, although the correspondence between the Vmpt" 
dependences of the DOS and T P is not perfect, probably due to 
the opening of the hybridization gap in the periodic Anderson 
model. 

As a consequence of the reduction in the f spin conduc- 
tance in P alignment, the TMR ratio dramatically reduces as 



FIG. 3: TMR ratio as a function of V imp f for e F = -1.2. Solid and 
open circles are the results for c imp = 0.1 and 0.2, respectively. 



Vmpt approaches e F , as shown in Fig. 3, where the results for 
Cimp = 0.2 as well as those for c; mp = 0. 1 are presented. The 
results indicate that the TMR ratio decreases when VBS ap- 
pears near the Fermi level and weakens the s-electron tunnel- 
ing. The reduction of the TMR ratio due to magnetic impuri- 
ties, however, is quantitatively stronger than the experimental 
value. The might be due to the rather strong s — d mixing 
chosen in this model. 

The relative position of VBS of magnetic impurities in tran- 
sition metals has been intensively studied before; 20-22 VBS for 
Cr impurities is close to the Fermi level, whereas that of V im- 
purities is higher than the Fermi level. For Cu impurities, on 
the other hand, V m pt may be nearly the same with = -2.5, 
and does not form any VBS. Therefore, we may conclude that 
impurities corresponding to Cr and V may bring about dra- 
matic reduction in the TMR ratio whereas Cu impurities do 
not. The results of both the conductance and TMR ratio are 
thus in good agreement with the experimental observations. 14 
Further, the present model showed no reduction in the con- 
ductance of the t spin state when magnetic impurities were 
placed on the Co leads several layers apart from the interface. 
Rather, it slightly increased near V; mp t ~ 0. This result is not 
inconsistent with the experimental result. 14 



IV. COMPARISON WITH CIP-GMR 

In order to compare the results obtained above with similar 
effects reported for GMR, we performed simple calculations 
of GMR in trilayers doped with magnetic impurities near the 
interface. As a simple model of Co/Cu/Co type trilayers, we 
adopted the s - d model, the Hamiltonian of which is given by 
adding the following Hamiltonian to eq. (2): 

H^-fJ^dldjcr, (7) 

ijcr 



because J-electrons may carry the current in this case. The 
n.n. hopping integral for af-electrons is taken to be 0.1. Now 
the (/-band is broadened by the <f-electron hopping, and we 
take a smaller value of the s — d mixing; y = 0.3. 

The lattice structure is assumed to be simple cubic, and the 
thickness of each layer is /6/12/6/ atomic layers in units of 
the lattice constant. Instead of the site representation used 
for TMR case, we use a mixed representation (/, ky) to treat 
the semi-infinite systems where I is an index of the atomic 
layers and ky = (k x , k y ) is a Fourier transform of the sites on 
each atomic layer. The Green's function G a (z) = (z - H a y l 
with z = £ + ir] is expressed by a tridiagonal matrix where the 
diagonal elements consist of 2 x 2 matrices of 

and the off-diagonal elements are alos given by 2 x 2 matrices 
including the hopping integrals of s- and <f -electrons. The d- 
level depends on the layer, and is expressed as \fi in eq. 
(8). As in the TMR case, we take uf T = -2.5 and ^ = -2.0 
when I belongs to the ferromagnetic (Co) layers, and = 
iff^ = -2.5 when I belongs to the nonmagnetic (Cu) layers. 
Because a set of parameter values different from those for the 
TMR case is used for the GMR case, the position of ep is 
different from that taken for the TMR calculation. 

To reproduce the GMR effect in our simple model, we in- 
troduce the interfacial spin-dependent scattering caused by an 
intermixing layer at interfaces. 23-27 Because the c/-level of the 
minority spin state of Co atoms is above the Cu J-band, Co 
atoms intermixed with a Cu layer will form VBS in the mi- 
nority spin state of the Cu s-band, resulting in spin-dependent 
scattering to produce the GMR. Here, we introduce the fol- 
lowing self-energy into the ^-component of the Green's func- 
tion: 



which represents the effects of Co atoms dissolved into the Cu 
layer. Here c indicates the concentration of such Co atoms and 
A is the broadening of cZ-levels due to the s — d mixing, which 
is taken to be A = 0.1, keeping in mind the width of VBS 
calculated in the first principles. 22 We further take rj = 0.001 
to reproduce the spin-independent resistivity. Here, we note 
that the concentration c itself does not have sufficient physical 
meaning, but cy 2 controls the magnitude of the self-energy. 
The conductivity parallel to the layers can be calculated using 
the Kubo formula without vertex corrections, because there 
is translational invariance along the planes. The GMR ra- 
tios thus calculated increase with increasing c and reach about 
60% fore = 0.5. 

After successfully reproducing the GMR effect in our sim- 
ple model, we introduce additional magnetic impurities into 
the adjacent layers of the intermixed interface. For simplicity, 
we use the same expression as that in eq. (9) for the self- 
energy for the majority-spin state of the s-band caused by the 
additional magnetic impurities, while c and iv are replaced 
with c; mp and Vi mpix , respectively. As in the calculation of 




FIG. 4: Calculated results of the GMR ratio as a function of the 
impurity potential V imp1 - for 6/12/6/ and 6/6/6 trilayers, respectively. 



TMR, we take Vi mp j, = and treat V; mp f as a variable pa- 
rameter. The GMR ratios calculated for c = c; mp = 0.5 with 
ep = -2.0 are shown in Fig. 4 (open circles) as a function 
of the impurity potential Vi mp f . It can be seen that an inverse 
GMR can be realized when V; mp f ~ 6f. Open squares in- 
dicate the results for a trilayer with 6/6/6 atomic layers with 
6f = —1.9, and the results show that the qualitative features 
are not affected by the thickness of the trilayer. The GMR ra- 
tio becomes minimal when Vimp-f ~ which characteristic is 
also found in TMR systems. The reduction, however, may be 
brought about by an increased scattering of s-electrons in the 
majority spin band due to s-d mixing with the impurity states 
in addition to a reduction in the s -component of the DOS. The 
effects of magnetic impurities such as Cr and V in the GMR 
case are almost the same as those obtained in the TMR case, 
though the origin of the reduction in TMR ratio is due to the 
blockade of s-electron tunneling caused by the anti-resonance 
in the majority spin state of s-band. 



V. CONCLUSION 

In conclusion, by a qualitative and systematic study of the 
dependence of the MR ratio on the type of impurities at the 
interface, we showed that both TMR and GMR ratios are af- 
fected strongly when the magnetic impurities form VBS in 
the majority spin state near the Fermi level. The results are in 
good agreement with the observed ones. Although the effects 
of magnetic impurities such as Cr and V are quite similar on 
TMR and GRM, the reduction in TMR ratios is due to a block- 
ade of i-electron tunneling, whereas that in GMR is mainly 
due to an increase of s -electron scattering in the majority spin 
state. More realistic and self-consistent calculations, however, 
are desired to obtain a quantitative understanding between the 
MR ratio and the interfacial electronic states. 
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